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Abstract 

The Eulcr-Gauss linear transformation formula for the hypergeo- 
in ■ metric function was extended by Goursat for the case of logarithmic 

singularities. By replacing the perturbed Bessel differential equation 
by a monodromic functional equation, and studying this equation sep- 
•^r , arately from the differential equation by an appropriate Laplace-Borel 

r) • technique, we associate with the latter equation another monodromic 

relation in the dual complex plane. This enables us to prove a dual- 
ity theorem and to extend Goursat 's formula to much larger classes of 



CN 



■S, 

^ . functions. 

1 Introduction 
> 

f^ . We propose a new approach for the study of the perturbed Bessel differential 

jn I equation (pBde). This approach is based on the duality between the func- 

tional monodromic relation generated by the pBde and its Laplace-Borel 
^ dual. The solutions of the latter relation can be viewed as generalized hy- 

^^ , pergeometric functions in the sense that they inherit a number of important 

CN '_ properties of the classical hypergeometric functions. 

The pBde, with an infinite number of parameters, gives rise to a mon- 
odromic relation with only two parameters 



P (Ce"') = P {Ce'^') + Te-<P (C) , (1) 

where a > 0, T is a complex constant which depends on all parameters of the 
pBde, P (C) is analytic at every point C 7^ and bounded at infinity in any 
sectorial region. The dual function F (i) is analytic in the t-plane cut along 
the interval (— oo, — a) of the real line and admits an analytical continuation 
along any path not containing points t = and t = —a. Among other 
properties which follow from ([T]) are that the function F (t) has exponential 
growth at infinity of minimal type in any sectorial region, and satisfies the 
monodromic relation 

F {{t + a) e^^ -a) =F{{t + a) e'^' -a) -TF (- (t + a) e^*) . (2) 



The essence of our approach is in considering both the monodromic re- 
lations for the pBde and its dual quite separately from the differential equa- 
tions which generate them. In brief, given a and T, we consider the linear 
spaces S and H of all functions retaining the analytic properties stated in 
the previous paragraph, and satisfying ([1]) and ([2]), respectively, and prove 
the following result: 

Duality Theorem. For fixed values of a and T, the Laplace transform 
operator C : H ^ S is a bijection. 

The monodromic relation ([5]) does not, on its own, determine the be- 
havior of F (t) at i = —a. However, the relation ([1]), together with the 
conditions on P (Q stated above, shows that F (t) has a logarithmic singu- 
larity at t = —a. Moreover we are able to prove the stronger statement: 

Theorem 2. For fixed values of a and T, assume that F (t) belongs 
to the linear space H. Then in the region |t -|- a| < a, |arg (t + a)\ < vr, the 
following relationship is valid 

oo 

F{t) = -TF{-{t + a))log{t + a) + Y.^k{t + a)\ (3) 

fc=0 

where the series in ^ is absolutely convergent and A/^ are complex coeffi- 
cients which can be evaluated explicitly. 

Formula ([3]) can be viewed as an extension of the special case of the Gour- 
sat linear transformation formula for the hyper geometric function 2F1 (a, b; c; t) , 
valid when a -|- b = c = 1 and 2F1 (a, b; c; t) has a logarithmic singularity at 
t = 1. While Goursat's approach relied on properties specific to the hy- 
pergeometric function, we establish a linear transformation formula for all 
elements of space H. 

The above results allow us, for fixed values of a and T, to provide all 
functions from the linear space S with asymptotic expansions in the C-plane, 
and using appropriate exponentially small terms to supply these expansions 
with uniform error bounds, see [5J. Returning to the original differential 
equation and developing our technique further, we derive an explicit formula 
for the connection coefficient T, see [B]. 

2 Two Relations for the Hypergeometric Function 

In this section we consider the linear transformation formula and mon- 
odromic relation for the classical hypergeometric function 2-^1 (i) ^; c; t) . 
The Euler linear transformation formula which determines the behavior of 
the hypergeometric function 2-^1 (Qj ^; c; t) at point t = 1 is well known (see 
for example pjj, 15.3.6) and is valid for all values of the parameters except 



for cases a+ b = c + m, where m is an integer. The associated formula 

2F1 (a, b; c; 1 - (1 - t) e^^-) - (2F1 (a, b; c; t)) 
= r±(a, b, c) (1 - t)^-"-" 2i^i (c - 0, c - b; c - - b + 1; 1 - t) , (4) 

where 

is less known and can be derived from the Euler linear transformation for- 
mula (see for example [7], section (4.2.4) formula (25). We refer to this 
formula as the monodromic relation for 2-^1 (a, b;c;t) . It is interesting to 
note that formula dU is valid for all values of the parameters including cases 
a + b = c + m referred to above. For the particular case a + b = c, this 
formula takes the form 

2F1 (0, b; + b; 1 - (1 - t) e^^^^) - 2i^i (a, b; a + b; t) 

= ^ ^^ v7r)r(lV ^^ ^''' ^' ^' ^ " *) ■ (^) 

In 1881, E. Goursat, [8], [9], when studying the logarithmic singularities 
of the hypergeometric function 2-^1 (d; b; c; t), derived supplementary linear 
transformation formulas for the cases a + b = c + m where m is an integer, 
which are quite different from Euler formula given in [Tj, 15.3.6. In partic- 
ular, he proved that if a+ b = c then for | arg(l — 1)| < vr and |1 — 1| < 1 the 
following formula is valid, ([Ij, 15.3.10,) 

2Fi(a,b;a + b;t) 

= i^^ (f ^7^<'(") (1 - ')" - '°«5 (1 - *) .F, („, .; 1; 1 - t)) , 

(7) 

where 

d{n) = 2^ (n + 1) - V (a + n) - V (b + n) , -0 = (log T)' . (8) 

Note that formula ([7D also yields the monodromic relation 1^. But the task 
of deriving ([7]) using only the monodromic relation ([6]) seems difficult, if not 
insurmountable, even for the particular case when a + b = 1. In this latter 
case formula ([6]) takes the simpler form 

2F1 (a, 1 - 0; 1; 1 - (1 - t) e±2-) - 2^1 (a, 1 - 0; 1; i) 
=r(o)2Fi(o,l-a;l;l-t), (9) 



where 

T (a) = =F2i sin (vra) , (10) 

and the same function occurs in each term of the relation. Even for this 
special case formula ([7]) remains as complicated as before. 

Given a > 0, replacing t — )• — - in ([UJ, this relation for the hypergeo- 
metric function 2-^1 (o, 1 — 0; 1; — ^) can be rewritten in the form ([2]), where 
T is given by pU|) . In the next section we show that the dual equation ([T]) 
for ([2]) arises from the Bessel differential equation. It turns out that the 
monodromic relation generated by the pBde is indistinguishable from the 
analogous relation generated by the Bessel differential equation. Therefore 
by identifying additional minimal properties of 2-^1 (ci;^;c;0 that are re- 
quired to derive ([7]) from ([6]) for the case a + b = c =1, or Bessel's case, we 
find similar conditions for the monodromic relation for the pBde. 

While our Theorem 2 appears to involve the derivation of d?]) from Q, 
assuming a+ b = 1, it in fact provides a partial answer to a deeper question 
regarding the additional properties of 2-^1 (ci, b;c;t) required to derive ([7]) 
from ([6]) . To answer this in full the results presented here for the pBde would 
need to be extended to a more general case of the perturbed Whittaker 
differential equation (pWde). 

3 The Spaces Sa^r and i/a,r 

Given a > and a complex number T, we define two function spaces Sa,T 
and Ha^T in terms of the monodromic relations ([1]) and ([2]) ■ 

3.1 Monodromic Relations generated by the pBde 

We begin with a differential equation of the form 
(Pu (a} oo a\ 



where a > 0, ao,ai, . . . are complex numbers, and the series is convergent 
for any complex C 7^ 0- Thus, if we set 

A{C) = ao + ^ + ..., (12) 

then A (C~^) is an entire even function. Equation (jlip can be considered as 
a perturbation of the standard Bessel equation. Indeed, upon applying the 
transformation ^ = iz,u = ziy^ Bessel's equation 

y" + V+(l-^)y = (13) 



reduces to the following particular case of pT]) with a = 2 and A{() = v 
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Clearly, ([TT]) is invariant under the rotations C — )• (e^'^^ and also un- 
der the reflection C ^ ~C- Since C = and C = oo are the only singular 
points of (jlip , any solution of this differential equation admits an analytical 
continuation in the C-plane along any path not containing these points. 

Let u (C) be a solution of ([TT]) decaying on the positive ray. Representing 
this solution in the form 

^(C) = e-Kp(C), (15) 

it is well known that the function P (Q has the same limit po as ^ ~^ oo 
along any ray in the region — ^ < arg^ < ^. 

Any pair of solutions consisting of u ((") together with one of the functions 
u (— C) , u (Ce'^*) and u (Ce~'^*) forms a basis in the space of solutions of ([TT]) . 
So, for example, we have 

u (Ce-*) = Au (Ce--^) + Bu (C) , (16) 

where A and B are constants. Using (J15p we can rewrite the relation ()16p 
in the form 

P (Ce^^) = AP {Ce-^') + Be-^'^P (C) . (17) 

Noting that all the functions P (Ce~'^') , P (C) aiid P (Ce'^*) tend to po as 
C -^ +00 and letting C -^ +oo in ([TT]), it follows that A = 1. Setting B = T, 
we rewrite the relation (|17p in the following final form 

P (^e^^) = P (Ce'^0 + ^e-'^'^P (C) , (18) 

where 

r = T(a,ao,ai,...) (19) 

is a complex constant which is usually referred to as the connection coeffi- 
cient, or Stokes multiplier. For the particular case of (jlip given by (|14p the 
constant T is known, 

T = 2i cos z^vr. (20) 

In the class of all differential equations with an irregular singular point 
at infinity of Poincare rank 1 the pBde should be regarded as an exceptional 
or degenerate case in the sense that a pBde generates a pair of monodromic 
relations which are not intertwined (each relation involves only a single func- 
tion) and the second relation can be derived from the first. For example, 
the relation for P {—Q follows immediately from (|18p . 

Representing P {() as a Laplace transform 

P{C)=C{F}:=C e-^'Fit)dt, (21) 

Jo 



it can be shown using the standard Laplace-Borel technique that F (t) sat- 
isfies the nionodromic relation 

F {{t + a) e^^ -a) =F{{t + a) e"" -a) -TF (- (t + a) e^^) . (22) 



The relation (j22|l can be viewed as dual to (|T8|1 . Moreover, it can be proved 
that F (t) has all the other properties listed in the second paragraph of 
Section 1 and also satisfies relation Q. In particular, for Bessel's case given 
by ()14p we have 

F{t)=2Fi(^^-u,^ + u,l;-iy (23) 

and as we noted in Section 2 the above properties are well known for this 
function. Analysis of the proof of the validity of these properties shows that 
instead of the pBde, we actually use the functional equation ([TH]) and the 
above properties of P (C) which have been extracted from (|11|) . 

Considering the relation (jlSp separately from the differential equation 
(|11|) . it turns out that (|18p . while being very much simpler than (|lip . yet 
remains a very rich source of information. 

3.2 The linear space Sa,T 

Given a > and complex T we introduce a linear space Sa,T of functions 
P (C) satisfying the conditions: 

i. P ((") is analytic in the ("-pla-ns punctured at the points (" = and 

C = oo; 

ii. Given P (C) G Sa,T, there exists a positive decreasing function Mp (r) 
such that if |argC| < tt and |(^| > r for any r, < r < oo, then 

\P {C)\ < Mp (r) . (24) 

iii. Every P {() G Sa,T, is a solution of the equation (fT8]l . 

Definition 1. Given P (C) G Sa^T we call the relation il8\) the Stokes 
monodromic relation. 

If P (C) is given by ([II]) and ([l5|) and T is given by ^ then P (C) G S'a,T- 
This shows that for any pair (a, T) there exists a non-trivial P ((^) G Sa,T- 
For instance, given T with u defined by (|20p . we represent a solution u [C) = 
Uy (C) of the Bessel equation ([HD decaying on the positive ray in the form 
u, (C) = e-^P, (C) . Clearly P, (f C) G Sa,T- 

Despite the fact that Definition 1 was suggested by a differential equa- 
tion, in order to distinguish those techniques which depend on the equa- 
tion itself from those which do not we have introduced a more general 
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object. The space Sa^T with independent constants T and a may con- 
tain elements not associated with such a differential equation. Indeed, it 
is easy to check that given an entire function $ (z) and P (C) E Sq r we have 
^{l/e)P{C)eSa,T- 

It follows from the definition of Sa^T that any element P ((") G 5a,r 
remains bounded at infinity in the region |arg("| < ^, and that in any 
sectorial region of the C-plane the function P {(,) is of exponential growth 
at infinity with exponent a. More precisely, we claim that the following 
statements are valid. 

Starting with Si = 1, Ti = T we introduce two sequences {5^}^^ and 
{^fclfc^i which are determined by the following recurrence relations 

S2m+1 = S2m, 5'2m+2 = S2m+1 + T2m+lT, (25) 

T2m+1 = T2m + S2mT, T2m+2 = T2m+l- (26) 

Lemma 1. Let P (C) G Sa^r o-nd let sequences {Sk}'^^i and {Tk}'^^i 
he defined by Ii25\) and i26\). Then for all integers m > —1 the following 
monodromic relations are valid 

P ((6(2-+!)-) = S2m+iP (Ce-"') + T2„+ie-<P (C) (27) 

and 

P (Ce(2™+2)^A ^ g^^^^p (^) ^ T2m+2e'"^ {P {(e~^')) ■ (28) 



These relations can be extended to all negative values of k. 

Illustration. One can check that for the Bessel case ()14p the relations 
(f27|) and ([28|) with coefficients satisfying ([25]) and (f26j) are identical to the 
formulae for analytical continuations of Hankel's functions (4.13) and (4.14) 
from [2], respectively. 

Proof. We will use mathematical induction on m to prove that the re- 
lations ([27|l and (f28|) are satisfied for some coefficients S2m+i, T2m+i, S2rn+2 
and T2m+2 , and in doing so we obtain the recurrence relations ([27]) and ([28]) . 

For 7TT- = — 1 we have the trivial instance of ()27p . with ^-i = 1 and 
T_i = 0; and the trivial instance of (j28p . with Sq = 1 and Tq = 0. For 
m = ()27p follows immediately from ()18p . along with the initial conditions 
5i = 1, Ti = T. From (fTSj) we also have, using (fTHj) . 



p (^e^-^) = P ((Ce"^) e"^) = P ((Ce"^) e-"*)+re-''^""> (Ce"*) = P (O+Te^^P (Ce"^) 

so that 

p [(e'^^i) = P (C) + Te'"^P (Ce'^^) . 

Now applying ([18]) again gives 



p (^e^^^) = P (C) + Te"^ (P (Ce"''^) + Te'^'^P (C)) 



which yields finally 

p [Qe^^i) = (1 + T^) P (C) + Te<P (Ce"'^*) . 

This shows, incidentally, that 5*2 = 1 + T^ and T2 = T, and so the recurrence 
relations hold when m = 0. 

Suppose now that odd and even hold for some value of m; we can use a 
similar process to prove that they also hold for m + 1. Thus, using ()28p we 
have 

p (^^g(2(™+i)+i)^*^ ^ s^^^^p (^g^i) + r2™+2e'^^^"> (C) 

Applying (fT8|) gives 

p (^^^(2m+^).^^^ ^ g^^^^ f^p (^g-.i) ^ ^g-aCp (^)^ ^ T2^+2e-'^^P (C) 

from which we obtain 

p f^^(2m+3).A ^ g^^^^p (^^-.^^ ^ (j.^^^^ ^ 52„+2T) e-<P (C) , 



and also the corresponding recurrence relations. 
In a similar way it can be shown that 

and thus that the last pair of recurrence relations hold. ▲ 
The next statement is an immediate corollary of Lemma 1. 
Lemma 2. Let P (C) £ Sa,T- Then for every sector 

S{r,e):={C:\C\>r,\aTgC\<e} 

there exists a positive constant Mp {6, r, T) such that for (^ G S {r, 6) 

|P(C)l<Mp(0,r,r)exp(a|C|). (29) 

Proof. The validity of 1^ follows from relations dS?]), (I2SD and the in- 
equality dHD.A 

We note that every P {Q G Sa,T can be provided with an expansion 
of the form 'Yl'k=oPkK^ ^ with complex coefficients pk, valid in the region 
— ^ < arg C, < ^. Introducing the remainders 

n~l 

P„(C):=P(C)-^S'" = 0'1'---' (30) 

fc=0 ^ 

it is possible to prove that 

Theorem 1. For P (Q € Sa,T the following estimates are valid for every 
n = 0, 1, . . ., and for \C\ > r 



TT 77-' 

largCI < - ^ |P„(C)I < Mp{r) -^^ (31) 



largCI < vr ^ |P„ (C)| < Mp (r) ^'^^ , (32) 



-<l-gCI<f ^|Pn(C)l<2Mp(r)^:|j|±^, (33) 

where Mp (r) is a constant which is proportional to that from ([2 
It follows immediately from ()32p that for n = 0, 1, . . . 



Tl 

\Pn\<Mp — ,Mp= inf Mp (r) . (34) 

a" 0<r<oo 

The proof of Theorem 1 is based on the techniques of the paper [3], 
Section 5, using the integral representations (110) and (111) given there and 
Theorem 2 of the current paper. We will not prove this in detail here since 
a more general result will be presented in our paper [5j. 

The reason for the deterioration when passing from (j32p to (j33p is the ap- 
pearance of exponentially small terms upon crossing the Stokes rays arg (^ = 
zbvr. This, which is a manifestation of the Stokes Phenomenon, will be dis- 
cussed briefly in the next section. 

We note finally that Lemma 1 and Lemma 2 allow us to extend the error 
bounds for P {() given by ([52]) to the whole Riemann surface — oo < arg ( < 
oo by adding to the remainders appropriate exponentially small terms upon 
crossing the Stokes rays, argC = zizrmr, m G N. 

3.3 Exponentially small terms and Stokes' Phenomenon 

The relation (jlSp provides an immediate clarification of the role of exponen- 
tially small terms in Stokes phenomenon, as we will now show. 

Comparison of relations (|32p and (|33p shows that upon crossing the rays 
arg C = zLtt the approximation given by ()32p begins to deteriorate and the 
greater the deviation from the ray the greater this deterioration becomes. 

The cause of the above deterioration is hidden in relation (jlSp . Indeed, 
relation (|18p generates the system of relations 

P„ (Ce-^) = P„ (Ce-^O + Te-<P {() , (35) 

where n = 0, 1, . . . and Pq (C) = P (0- L^t us introduce functions 

E+ (C) = Te<P {(e-^') , vr < arg C < y (36) 

E- (C) = -Te<P (Ce^O , -^ < argC < -vr. (37) 



These functions decay exponentially in the regions { C : ^r ^ arg C < t" } ^^'^ 
{C : — ^ < argC < — vr} , respectively, as (" — )• oo. Using (j35|) we observe 
that subtracting these exponentially small terms from the remainders, we 
retrieve the accuracy of (p2]) : 



^ < argC < ^ ^ |Pn (C) - E+ (01 < Mp (r) ^^^^, (38) 

-f < argC < -vr ^ |P„ (0 - i^- (C)| < Mp{r)^^^^, (39) 



Indeed, it follows from ([35]) . ([36]) and ([37|) that for n = 0, 1, . . . 

7r< arg C < ^ ^ arg (Ce'^-) G [-^, tt] : P„ (C) - E+ {() = P„ (Ce'^-) 
-< argC < -^ ^ arg (Ce^-) G [-^,7r] : P„ (C) - E" {() = P„ {(e^^') , 



which, using ([32]) yields (p8]l and (p9]) . 

Note that if the monodromic relation is generated by the differential 
equation (jlip then evaluation of the exponentially small terms reduces to 
evaluation of the constant T. To clarify the cause for the above phenomenon 
we must turn to the dual complex plane and to consider a dual space Ha^T 
of functions F (t) analytic in the dual plane. 

3.4 The linear space Ha^x 

Given a > and complex T we introduce a second linear space Ha^T of 
functions F (t) satisfying the conditions: 

• (i) F (t) is an analytic function in the complex t-plane punctured at 
the two finite points t = and t = —a, and with exponential growth 
of minimal type at oo in any sectorial region {t : — oo < a < arg t < 

/3 < +oo}; 

• (ii) there exists a branch of F{t) which is analytic in the t-plane cut 
along (— oo, —a] ; 

• (iii) the branch given by (ii) satisfies the estimate F(t) = O (log (t + a)) 
as t — )• —a, |t + a| < a, |arg (t + a)| < vr; 

• (iv) the branch given by (ii) satisfies the monodromic relation (|22p 
given by 

F ((t + a) e^* -a)=F{{t + a) e'^' -a)-TF (- (t + a) e^^) , (40) 

where —a < t < +oo , and T is the constant from (fTSl) . 
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Note that the points (t + a) e'^* — a and (t + a) e~'^* — a belong to the 
upper and lower banks of the cut, respectively, and the relation (j^Oj) 
can be rewritten as 

F (te^* -a) =F (ie"" -a) -TF (-te^^) , i > 0, 

as well as in the equivalent form 

F {{t + a) e-2-^ -a)=F{t) + TF (- {t + a)) , (41) 

where t belongs to the upper bank of the cut. 

Definition 2. Given F (t) S Ha^T we call the relation ^JW the hyperge- 



ometric monodromic relation. 

Remark 1. Since function F (t) is multi-valued, it is necessary to inter- 
pret the expression F (— t) "carefully". If Fq (t) is a branch of F (t) which is 
given by (ii) of Definition 2, then Fq (— t) is a single-valued analytic function 
in the t-plane cut along [a, +oo) such that 

\Fo(te-^*) if9t>0. 

We define F (— t) as an analytical continuation of Fq (— t) to the complex 
t-plane punctured at the two points t = and t = —a. The function 
F {—t — a) satisfies condition (i). 

Below we will also use F (t) to denote the branch given by condition (ii), 
where this will cause no confusion. 

For the case given by ([23]) setting F^, (t) = 2-^1 {I ~ ^^l + ^A',-^), 
where a = ^ — u and T = 2i cos utt, one can check that F^, (t) G Ha^T- 

In the next section we will show that there is a one-to-one correspondence 
between the linear spaces Ha^r and Sa,T- 

4 The Duality theorem 

Defining the Laplace transform operator C by (|2ip we claim that 

Theorem 2. Given a and T, the operator C : Ha.T — ^ Sa,T is a bijection. 
We present the proof in two parts. 

4.1 Part I: CHa,T C Sa,T 

Proof. Let F (t) £ Ha,T and let P (C) be given by ([21]). Then it follows 
from conditions (i) and (ii) of section 2.3 that P {C,) admits an analytical 
continuation from the right half-plane to the C-plane cut along the interval 
(— cxD, 0) as a function bounded at 00 and continuous in the cut plane except 
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for C = 0- Moreover, for — vr < 9 < tt and for < p < +c«, using Cauchy's 
theorem, the integral representation 

P (^pe'^^ = pe'^ / e-P^^'^F {t) dt (42) 

is vahd, where the integral is absolutely convergent. Setting 6 = ibvr in 
([12|) . we retain the previous convergence due to the condition (iii) of 2.3. 
Subtracting the second integral from the first we have 

P {pe""') - P {pe-^') =-p eP^F {t) dt + p e^*F (t) dt. (43) 

Jo Jo 

Due to condition (ii) F (t) is analytic in the circle \t\ < a, and thus for 
t G (-a, 0) we have F {t) = F (t*), where t* = {t + a) e'^'"* - a. Therefore, 
writing J_ for an integral along the upper bank of the cut, ()43p can be 
rewritten in the form 

/~oo /•— oo 

eP^F {t*)dt + p eP^F (t) dt (44) 

-a J ~a 



-oo 



= -p eP'{F{f)-F{t))dt. (45) 

J —a 

Using the monodromic relation ()40p in the form ()4ip . the last integral in 
()45p can be rewritten as 

eP^F (- (t + a)) dt = Tpe^^P / eP'-^+'^^F (- (t + a)) dt. (46) 

-a J —oo 

Finally we have 

/O r+OD 

eP^F (-t) dt = Tpe-^P / e^^^'F (t) dt = Te-^PP (p) . (47) 
-oo Jo 

The chain of relations ()45p - ()47p shows that P (C) satisfies the condition (iii) 
of definition 1, and thus, P {() G Sa^r-^ 

4.2 II. Sa,T C CHa,T 

We assume now that P (C) G Sa^T- We set -F (t) equal to the Borel transform 
of P (C) , so that P (C) = C {F} , and prove, in five steps, that F (t) G Ha^r- 
(1) Integral representation for F (t) for |argt| <^. Given r > and 
^ > 0, we introduce the contour 751 (r) with an anti-clockwise orientation as 
follows: 

le (r) = l-e (r) U Ce (r) U h (r) , (48) 

where 

leir) = {C:C = pe'',r<p<oo] (49) 
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and 

Ce{r) = {C:-9<a.TgC<e,\C\=r}. (50) 

For t > the function F (t) can then be represented in the form 

F(i) = ^l «*^^(C)^, (51) 

2vri J^^(r) C 

with any f < 6* < ^ and r > 0. 

Below we set 9 = n and omit the subscript vr in Ctt (r) and in 7,^ (r) if 
this does not lead to confusion. 

First we prove that this integral converges for t > 0. Then we demon- 
strate how, by changing the integral representation, we can obtain an ana- 
lytical continuation of F (t) that allows us to check the validity of conditions 
(i)-(iii) of 2.3. Finally we show that the function F (t) satisfies the relation 
(liOD . Setting 

(52) 
(53) 



and 



we have 



F+it,r) = 


2vrz Ji^^r) C 


F-{t,r) = 


27ri Jl_^{r) C 


Fo{t,r) = 


2vrz Jc{r) C 



(54) 



F (t) = F+ (t, r) - F- {t, r) + Fq (t, r) . (55) 



Observe, that Fq (t, r) admits an analytical continuation from the positive 
ray to the whole i-plane as an entire function of t. Since P (C) G Sa,T, the 
relation (|24|) yields the following estimate valid in the whole t-plane 

\Fo{t,r)\<Mp{r)e^^^\ (56) 

meaning that Fq (t, r) is an entire function of exponential type r. On the 
other hand, for the functions F^ (t, r) given by (|52p and (j53p the inequality 
(|24p yields the estimates 

|f± (,,)| < ^ /*~e-C-,, = ^,-,0 < * < +00, (57) 

The estimates ([56]) and ([57|l show that integrals (f52]l - ^^ are absolutely 
convergent, so that the function F (t) is well defined for positive t, and that 
its growth on the positive ray is determined essentially by estimate (|56p . It 
follows from (|56p and the preceding comments that the integral in (j2ip is 
absolutely convergent for 3f? {(^} > r. Since r is an arbitrarily small positive 
number it follows that this integral is absolutely convergent for all C such 
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that |arg (\ < ^. Therefore formula ()5ip can be considered as the converse for 
the Laplace transform ()2ip , which proves the validity of the representation 

e-<F(s)ds,--<argC<-. (58) 

(2) Representation for F (t) in the region {t : \t\ > a} . Let 9 be 
any angle, — oo < 9 < oo. To obtain the desired analytical continuation of 
F (t) we rotate the path of integration in ()5ip . with 9 = n, about the origin 
through an angle —9 while simultaneously increasing argi by 9. Lemma 1 
yields 

|P(C)|<M(0,r,T)exp(a|C|), (59) 

where the positive parameter M {9,r,T) depends on 9, r and T. Following 
the above rotation the integral in the form 

^M^=^/" e*^"V(C)^,0<t<oo, (60) 

retains the absolute convergence of the integral (j5ip . Thus, using Cauchy's 
theorem, it can be shown that ()60p provides an analytical continuation of 
(j5ip . Since 9 is arbitrary, this shows that F (t) admits an analytical contin- 
uation to the exterior of the circle {i : |t| = a}, and ()60p provides an integral 
representation for F (t) in the region {t : |argf| < vr} and also in the region 
{t : |t| > a, — oo < argf < oo} . It also follows from ([5UP and ([^ that given 
< r < oo in any sectorial region {t : \t\ > R > a, a < argt < /3} , of the 
t-plane the following estimate is valid 

|F(t)| <i^e"l*l, (61) 

where the positive parameter K depends on R, a, (3, and r. 

(3) Representation of F (t) in the region (f : |arg(i + a)| < |). 
Next we prove that F (t) is analytic in the half-plane ?R-Z > —a and therefore 
it satisfies the conditions (i) and (ii) of 2.3. In order to prove this we first 
assume that < t < oo and derive another integral representation of F (t). 
Now, returning to the representation (j5ip and making, in (j52p and (|53p . the 
change of variable suggested by (flQ]) for 9 = ibvr, it follows that (f5T|) can be 
written in the form 



F (t) = Fo (t, r) + — J e-*? (P (Ce") - P (Ce— )) ^, (62) 

where Fq {t, r) is given by (fM|) . 
Using the relation 

P (Ce*^) - P (Ce~*^) = Te~<P {() , 
14 



which fohows from ()18p , we can represent -F (t) , t > in the form 

rp r-+ao J/- 

F it) = Fo (t, + ^1 e-(*+")^P (C) -f. (63) 

Formula ()63p shows that F (t) can be continued analytically from the positive 
ray to the half-plane |arg (t + a)\ < ^. This completes the proof of step (3). 
Steps (l)-(3) prove that F (t) is analytic in the whole t-plane punctured 
at t = and t = —a, and formula (j63p represents a branch of F (t) which is 
analytic in the plane cut along the interval (— oo, —a]. We keep the notation 
F (t) for this branch, if it does not lead to confusion. Clearly the function 
F (t) is analytic in the circle \t\ < a. Applying Watson's lemma to the 
relation ([58]) we can provide every element P (C) G Sa,T with an asymptotic 
expansion of the form '}2T=oPk/C^^ where 

Pk = F^^^ {0),k = 0,l,...,F = C-'P. (64) 

Analytical properties of F (t) given by (l)-(3) show that the relation 

n-l .X 

^(C)-Ep^/C' = - ,C^oo, (65) 

A;=0 Vj / 

is valid for every n G N and for (^ satisfying — f < arg C < ^ , and can be 
extended to every sub-sector of the region — ^ < arg ( < ^. Next we use 
this fact to prove that F (t) satisfies condition (iii) of 2.3. 

(4) Behavior of F (t) at t = —a. We prove the following statement. 

Lemma 2. Let P (^) G Sa^T cind F (t) be given by II63\) . then there 
exists a constant Aq such that in the region \t + a\ < a, |arg (t + a)\ < it the 
following formula is valid 

F it) = -TT^ log {t + a) + Ao + oil),t^ -a. (66) 

ZlTl 

Proof. To prove ([66|) we use ([63]) to estimate F (t) . Note that Fq (t, r) 
is an entire function in the t-plane and from ()54p we may write 

Po{t,r) = ^^ I e('+-^<e-<P{0^. 

C(r) 

Since e(*+'')'^ = 1 + O (t + a) as t -^ -a it then follows that 

i^o(t,r) = ^ / e~<P{0'^+O{t + a),t^-a. (67) 

C{r) 

Let us denote the integral on the right hand side of (j63p by / (t, r) and 
observe that it can be written in the form 
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where po is given by ()64p with A; = 0. Using the relation (|65p with n = 1 
together with e^*"^")^ = 1 + O {t + a) as t —?• —a, the first integral in ([68]) 
can be represented as 



(69) 
where the integral on the right hand-side is absolutely convergent. 

Finally we evaluate the second integral in (j68p using the exponential 
integral 

El (z) = / e-^-— (70) 

Jl T 

and its power series expansion, see formula (5.1.11) in [1], 



oo 






i?i(z) = -logz-7-E r \ , V (71) 

^-^ (m + 1 ml 

m=l 



where |arg2:| < vr and 7 is Euler's constant. Thus, replacing z in ()70p by 
r {t + a), and applying ([7T]) we have, for |arg (t + a)| < tt, 

=^ f-log(r(t + a))-7- V ^-P^^^] . (72) 

27rz I ^^ ^ ^^ ' ^^ m + l)m! / ^ ^ 

\ m=l / 

Combining ([67|) . ([69]) and ([72]) . we obtain the following estimate for F (t) 

F{t) = -^r^log{t + a) + Ao + {t + a) ,t ^ -a,\aTg{t + a)\ < tt (73) 

2m 

where 

C(r) 

(74) 
and r G (0, 00) . While the individual terms on the righthand side of (|74p are 
clearly dependent on r it follows from ([73p . taking the limit as t — )• —a, that 
the coefficient ^0 does not depend on r. A 

It remains to be proven that F (t) satisfies the monodromic relation ([40p . 

(5) Derivation of the monodromic relation for F (t) . The proper- 
ties of F (t) given by (i) and (ii) of the section 2.3. which have been proven 
in the preceding paragraph allow us to use the relations ([^3]) and ([15]) and 
to show that for < p < +00 

/— 00 
gpi (F (t*) - F (t)) dt, 
-a 
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where t* = {t + a) e ^'^* — a. On the other hand using elementary transfor- 
mations given by ()46p and (j47p we have 

/— oo 
eP^F{-{t + a))dt. 
-a 

Since the P {pe''') - P {pe'''') = Te'P^P (p), it follows that F (t*) - F (t) = 
TF{-{t + a)). 

Steps (l)-(5) justify the inclusion Sa^T C CHa^x- A 

In the next section we use the Duality theorem to prove a version of the 
linear transformation formula for F (t) G Ha^T- 

5 Linear transformation formula 

The following statements is valid. 

Theorem 3. Assume that F {t) G Ha^r- Then for the branch of F {t) 
given by condition (ii) of 2.3 in the region \t + a\ < a, |arg (t + a)| < n the 
following formula is valid 

F{t) = J2 ^k {t + af - —F (- {t + a)) log (t + a) , (75) 

fc=0 

where the series in ( f75[ j is absolutely convergent, Aq is given by |7-^[ j, where 
Pq = F{0), and A^, k £ N, are complex coefficients which can be found 
explicitly. 

Remark 2. For BesseVs case given by ( |j^[ j relation ( (75^ follows im- 
mediately from a degenerate case of the Euler transformation formula for 
the hypergeometric function, see 15.3.10 of |I]/ and Appendix 1. Thus, the 
relation ((75| j can be viewed as a generalization for the elements of Ha.r of 
the linear transformation formula ^ with a + b = 1. 

Proof of Theorem 3. The relation ()66|) of Lemma 2 can be rewritten 
as 

F [t) = -^r^ log (t + a) + ^0 + o (1) , i ^ -a, (76) 

where ^o is given by (|74|) . Since F {t -\- a) is analytic in the circle |t + a| < a 
we have the relation 

F{-{t + a)) log {t + a) = {F (0) + O (t + a)) log (t + a) as t -^ -a, 

which yields 

F (0) log {t + a) = F{-{t + a)) log {t + a) + o (1) . 

Substituting the last relation into ()76p allows us to rewrite it as 

F (t) = ---F (- (t + a)) log {t + a)+Ao + o{l),t^ -a, (77) 
2m 
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Let us introduce a function ^ (t) such that 



T 
F (t) = -7—-^ (- {t + a)) log{t + a) + Ao + ^ {t + a) ,t ^ -a. (78) 

Clearly $ (t) is analytic in the t-plane punctured at the two points t = 
and t = a. On the other hand the relation (j77|) shows that 

$(t + a) = o(l),i^ -a. (79) 

Since |t + a| < a, |arg(t + a)| < vr we can continue analytically both 
sides of (|78p by rotating t + a through angles vr and — vr about the origin. 
This allows us to rewrite (ITHI) as 



F {{t + a) e'^' -a) =Ao + <^{{t + a) e^^) 
+^.F (- {t + a) e-*) log ((t + a) e^') , (80) 



and also as 



F {{t + a) e-^' -a)=AQ + ^{{t + a) e'^') 
-^.F (- (t + a) e-^') log ((t + a) e"") . (81) 

Subtracting ()8ip from (|80p and using the monodromic relation ()4ip gives 
F ((t + a) e^* - a) - F ((t + a) e"^* - a) = -TF (- (t + a) e^^) . 
A straightforward calculation using 

F (- (t + a) e^*) = F (- (t + a) e"^^) , |t + a| < a, 
shows that 
-TF (- (t + a) e^') = -^F (- [t + a) e"^) log {e^^')+^ {{t + a) e^')-^ {{t + a) e"-^) 

Canceling equal terms, the last relation can be simplified and we have 

$ {{t + a) e") =^{{t + a) e"") 

which, using ()79p . shows that ^ (t) is analytic and single- valued in the circle 
\t + a\ < a and that $ (0) = 0. Thus, 



$(t + a) = ^^fc(t + a)^^, 



fe 



$W (- 



■a 



k=l 
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6 The coefficients A^ 

In this section we describe a procedure that allows us to derive formulas 
for computing the coefficients Ak,k = 1, . . . , given by ([75]) . Simultaneously 
we give an alternative proof of Theorem 3 which is based on Theorem 1. 
Returning to the relation ()63p rewritten as 



F it) - Fo {t, r) = —j^ e-^t+aKp (^) A, (82) 

where Fq (t, r) is given by ()54p we represent the integral in (|82p in the form 

e-(*+a)Cp (^) A = j„ (t, r) + j„ (i, ^) (83) 



^^ f + OO 

2^1 



where 



and 






c 



(84) 



fc=0 -^ '^ ^ 

Given n G N, and expanding /„ (t,r), J„ (i,r), and Fq (t,r) in series of 
the form 

log {t + a)Y, Pk {t + af + Y^ ak {t + a)'' , 

k k 

we are only interested in the values of the coefficients a„ in each series. 
Let us denote these coefficients by a„ (/), a„ (J) and an (Fq) for /„ {t,r), 
Jn (t, r) and Fq (t, r) , respectively. Then ([82]) shows that the coefficient An 
at (t + a)" can be written as 



An = On (I) + a„ (J) + On (Fq) . 

Firstly, expanding the exponential e^^*^")'' in (|84p into a Taylor series 
we have 

(87) 
so that 



a„(/) = (-l)"^l 
2m n\ 



p+oo ( " \ 
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and the integral is absolutely convergent since [P {() — Ylk=oPk/C'') C" "^ = 
O i-pz) as C — 7- oo. Moreover, using error bound (j32p of Theorem 1 yields 
the estimate for n G N 



|a„(/)|<^^^^,Mp= inf Mp{r). (89) 

ZTT a" 0<r<oo 

Secondly, we calculate Jn {t, r) using the exponential integral 

En{z) = J e~^*-,n = l,2,..., (90) 

and its power series expansion 

(n — 1)1 ^-^ Im — n + ljml 

where 

n—\ ^ 

|argz| < 7r,V'(l) = -7,V'('^) = "7+ V" — , 

^-^ m 

7Tl=l 

and 7 is Euler's constant. 
We have 

'Jn{t,r) = —Y^E,^,{rit + a)) (92) 

k=0 



and substituting (|9T|) into (|92|) yields 

•^"(*''') = ^E^fe ^'^^t''^^ (-logr-log(t + a)+V'(fc + l)) 

/c=0 

2TTi 2^ r^ 2^ (m-k)m\ ' 
It follows from ([93]) that the coefficient a^ (J) for J„ (t, r) is given by 

aniJ) = —^-^Pn{-logr + i^{n + l)) - — ^^-^ -^. (94) 

2tti n! 2tti ^-^ n\[n — k) 

k=0 ^ ' 

Using ([M|) and setting r = 1, we have 

, , ,, TMp (1) V' (n + 1) , , ,, 

27r o" 

Finally, it follows from (j67p that 

1 1 /■ „-af„..^^C 



ar. 



n{Fo) = ^-, e-<P(C)^, (96) 



27ri n! ./cir) 
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and 

|a„(Fo)|< ^ (97) 

So we have proved that the coefficient An at (t + a)" is represented in the 
form given by (f86l) . where a„ (/) ,a„ (<^) , and a„ (Fq) are given by (f88l) . 
([M|) . and ([96]) . respectively. Analysis of the estimates for a„(/),an(J), 
and a„ (-Fq) given by (f89|) . (f95|) . and (f97j) shows that the term ([88]) gives the 
main contribution to the asymptotics of A„ as n — )• oo. We have the estimae 



, , , TMp(l) v/^r+3 , , ,, 

An = f^^ '^ 2 1 + 1 ,n^oo 98 

which shows that the power series X^n^o ^n^" i^ absolutely convergent in 
the circle of radius a. 

It turns out that we have an alternative proof of Theorem 3 based on 
Theorem 1. Indeed, combining expressions ([82]) . ([83]) . ([87|) . ([93]) . we have 
for |arg {t + a)\ < a 

fe=0 ■ ra=0 

(99) 
where ^„ are given by ([86|) . Since p^ = F*-^) (0) we have 

-ilf^^Eg ^ fc, ^M og(t + a) = --F(-(t + a))log(t + a), 

fc=0 

so (I99D can be rewritten as 



F (t) = -^—.F (- (t + a)) log (t + a)+y An (t + a)"+o ((t + af) ,t ^ -a. 

n=0 

It remains only to evaluate more accurately the remainder o((i + a)"') 
in ([87|) . Let us return to the expression (|8^ and represent it in the form 

C 



i.(.o^|.f^-^-)^(p(c)-|:wc^) 



which can be rewritten as 

»+oo 



^"^^''^^ = ^/ ~e-(*+"^^P„+i(C)^, (100) 

where P„+i (C) is given by ([HD]) . Using the expansion 

y.' (-ir_(t+arc 



e-(*+a)c = V (-ir^^+^rc- ^ (ziri+^rr ^-(,+.)c(o 



m=0 
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where r < C, (r) < c«, and notation ()88p . the expression (jlOOp can be rewrit- 
ten as 

n-l 

In{t,r) = Y,c^k (!) {t + a)'' + R{t,n,r) , 
fc=0 

where 

i2(t,n,r) = ^ ^"'^"l!'""^" /'"°°C"e-(^+")^(-)p„-,i (C) f • (101) 
Applying the estimate for the remainders given by (|32|) we have for n = 

|P„., (01 < MP (r) '" +4^ , 

which together with ()101|) yields the inequality 

,p., „ ^ TMp (r) (n + 1) ^^^+4 ,^ , ,„, 

i? t,n,r < \t + a\ ,n = 0,l,.... 

Using estimate (|86p this completes the alternative proof of Theorem 3. ▲ 

7 Appendix. The perturbed Whittaker equation 

We have replaced the differential equation (jlip by a system of functional 
equations 

PiiCe^') = Pi(Ce-"0+re-<Pi(C), 
P2(Ce'^0 = P2{Ce-''')+Te<P2{0. 

We note that these two equations of the system are not linked. A reason 
for this is that the coefficient A (C) is an even function. However, typically, 
such equations are intertwined. It is enough to add the term ^ to A (C) to 
obtain from (jlip an intertwined differential equation 

(Pu { a? h qq ai \ 

We consider a more general equation 

where B ((") = X^fc^o ^fc/C^ is an entire function of l/C- A special case 
a = 1,6 = — K, 6o = /i^ — 1/4, 6fc = for A; = 1,2,..., is known as Whit- 
taker's differential equation. Thus, ()102p can be considered as a perturbed 
Whittaker equation. Let us assume that a > 0. Then it can be proved that 



22 



there exists a pair of linearly independent solutions ui (C) and U2 (C) of ()102p 
such that 

ui (C) = e"t^C-«^i (C) , -f < arg C < '" 

n2(C)=etCc!p2(C),-f <argC<f 



2 --6S - 2 ' (103) 



where 

A(C),^2(C) = l + o(l) (104) 

as C — 7- oo along any ray of these sectorial regions. It can be shown that 
the solutions ui (C) and U2 (C) are uniquely determined by their asymptotics 
given by (|103p and (jl04p and that Pi (C) and P2 (C) admit analytical con- 
tinuation along any path not crossing C = 0- 

Since ui (Ce^'^*) and U2 (Ce^'^*) are also solutions of (|102p and since ui {(,) 
and U2 (Cc^'^*) are linearly independent solutions we have 



ui (Ce^-^) = Aui (C) + 5n2 {Ce^^') 
U2 (Ce2-^) = Cu2 (C) + Dui (C) 

where A,B,C,D are complex constants. Using (jlOSp . we have 



(105) 



n2(Ce2-^)=e2-^eKC"^2(Ce'"') 



which allows us to rewrite ()105p in the form 



e2-!et'^C^^2 (Ce^^*) = C7et^C^^2 (C) + Z^e-f^C"^ A (C) 

(106) 
Let us simplify (jl06p 

e-2"^«Pi (Ce^-*) = APi(C) + Pe2"^^e<C^i^2(Ce^"')' (10^) 
e2-^^P2 (Ce^-^) = CP2iC) + De-<C^Pi{0- (108) 

and analyze the last pair of equations. Considering the equation p07p we 
assume that argC e (-^, -f ) , thus arg (Ce^''*) G (f , ^) . Then letting C 
go to infinity, and taking into account that Pi (C) — ?• 1,-P2 (Ce^'^*) -^ l,and 

that the term Be"'^C~°- P2 (Cc^'^*) is exponentially small, it follows that 

A = e-^^'a. (109) 

A similar analysis for p08p and for arg(^ G (~f ' f) shows that 

C = e^^'a. (110) 
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Using (fT09]) and (fTTOl) . the system of equations (fTOTl) and ([TOH]) can be 
rewritten as 



Pi(Ce2-^) = P^iC) + Be^^'-ae<C^P2{Ce'^'), (111) 

P2(Ce^"0 = P2iC) + De-^^'-^e-<C^Pi{C). (112) 

Setting 



Ti =Be^'^'aT^ = De"27ri- 



we have finally 

Pi(Ce2-) = Pi(C) + Tie<C^P2(Ce'"0' (^3) 

P2(Ce'"0 = P2{Q) + T2e-<Q-"-tp,{0. (114) 

Note that the above pair of relations can also be rewritten as 

Pi(^e^i) = Pi(Ce-^^) +rie-2^^ae-<C^^2(Ce"), (115) 

P2(Ce"0 = P2(Ce~"')+^2e'"'«e"^r^i^i(Ce-"0- (116) 



Our principal idea is to consider the system of equations (jllSp and ()116p 
separately of differential equation. Assuming that Pi (C) and P2 (C) are 
analytic and bounded in the regions — ^ < arg C, <^ and — f < arg (^ < ^, 
respectively and representing these functions in the form 



Pi (C) = C / e-^*Pi (t) dt 
Jo 



and 

r>oo-e 



^oo-e 

P2 {C) = C e-^*P2 it) dt, 

Jo 



we claim the following statement. 

Theorem, (i) Functions Pi (t) and P2 (t) admit analytical continuation 
to the t-plane punctured at points t = and t = —a for Pi (t) and at points 
t = and t = a for P2 (t) ; 

(a) there exist branches of Pi (t) and P2 (t) analytic in the t-plane cut 
along (— c«, —a) and (a, +00) , respectively; 

(Hi) Pi (t) and P2 (t) satisfy a dual system of monodromic relations. 

The most nontrivial is the assertion (ii). Setting 

p«=p('=)(0),pl^)=pf(0),A: = 0,l..., 
it follows from this assertion, in particular, that 

-— < argC < — ^ lim Pi (C) = Pi (0) , 
-^ < argC < ^ ^ lim P2 (C) = P2 (0) , 
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and moreover that Pi {() and P2 {() can be expanded into asymptotic series 

Pk^ 37r ^ 37r 

^,-^<argC<- 

(2) t- 



00 (1) „ „ 

fc=0 

ft(o~i:C.-i<«gc<5f, 



fe=0 

8 Conclusion 

Our aim for the future is to extend the approach described here to more 
general systems of functional monodromic equations that are generated by 
linear differential equations or systems of differential equations with an ir- 
regular singular point of arbitrary Poincare rank at infinity. The pBde is 
the simplest differential equation that can be reduced to a single functional 
equation while preserving most of the difficulties arising in the general case. 
It is for this reason that in this initial study we have limited our attention 
to a detailed consideration of the pBde. 

9 Acknowledgment 

Authors wish to express their gratitude to Sergey Suslov for his useful sug- 
gestions. 

REFERENCES 

[1] M. Abramowitz, I. Stegun, (eds.), Handbook of mathematical functions, 
Dover Publ., Inc., New York, 1992. 

[2] F.W.J. Olver, Asymptotics and special functions. Academic Press New York 
San Francisco London, 1974. 

[3] V.P. GuRARii, V. Katsnelson, V. Matsaev and J. Steiner. How to 
use the Fourier transform in asymptotic analysis. Twentieth century harmonic 
analysis— a celebration (II Ciocco, 2000), (NATO Sci. Ser. II Math. Phys. 
Chem., 33), Kluwer Acad. Publ., Dordrecht, 2001, 387-401. 

[4] V.P. Gurarii. Error bounds, duality, and the stokes phenomenon, I Algebra 
and Analysis, 21, 6, 2009, 80-150, ( see, also AMS St. Petersburg Math. J. 21 
(2010), 903-956, S 1061-0022(2010)01125-7. 

[5] V.P. Gurarii and D.W.H. Gillam , Perturbed Bessel equations II. Error 
bounds and the Stokes phenomenon, (to be submitted). 



25 



[6] V.P. GuRARii, D.W.H. GiLLAM, AND V.I. Matsaev, Perturbed Bessel equa- 
tions III. Connection coefEcients. (to be submitted). 

[7] A. Kratzer and W. Franz, Transzendente Funktionen. Lepzig, 1960. 

[8] E. GOURSAT, Sur I'equation difFentielle lineaire qui adment pour integrate la 
serie hypergeometrique. Ann. Sci. Ecole. Norm. Sup. (2), 10, 3-142, 1881. 

[9] E. GouRSAT, Proprietes generales de I'equation d'Euler et de Gauss. Actualites 
scientifiques et industrielles 333, Paris, France, 1936. 

Mathematics 38, FEIS, 
Swinburne University of Technology 
PO Box 218 Hawthorn 3122; 
e-mail: vgurarii@swin.edu.au. 

School of Information Technology and Mathematical Sciences 
University of Ballarat 
P.O. Box 663 Ballarat, 3353, 
e-mail: dgillam@ballarat.edu.au 
VIC 
Australia 

VIC 
Australia. 



26 



